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Abstract 

A gas of particles which colhde inelastically if their impact velocity exceeds a certain 
value is investigated. In difference to common granular gases, cluster formation 
occurs only as a transient phenomenon. We calculate the decay of temperature 
due to inelastic collisions. In spite of the drastically reduced dissipation at low 
temperature the temperature surprisingly converges to zero. 
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1 Introduction 

Gases of dissipatively colliding particles (Granular Gases) in the absence of 
external forces reveal a variety of interesting phenomena, such as characteristic 
deviations from the Maxwell distribution [1-3], overpopulation of the high 
energy tail of the distribution function [4], anomalous diffusion [5,6], and 
others (see [7] for an overview). The most striking phenomenon which makes 
Granular Gases different from molecular gases is the self-organized formation 
of spatio-temporal structures such as clusters [8] and vortices [9] . 

The loss of mechanical energy of dissipatively colliding particles i and j is 
characterized by the (normal) coefficient of restitution, which relates the nor- 
mal component of the relative velocity before a collision, g, to that after the 
collision, g': 



with Vij = Vi — Vj and vlj being the corresponding after-coUision value. Fre- 
quently it is assumed that the coefficient of restitution is a material constant 
which simplifies analytical calculations considerably. This assumption, how- 
ever, contradicts experiments [10] and disagrees with a dimension analysis 
[11, 12]. Instead the coefficient of restitution is a function of the impact veloc- 
ity. 

Obviously, the coefficient of restitution can be obtained by integrating New- 
ton's equation of motion for the collision with appropriate choice of the contact 
force. The clastic component of the contact force of spheres of diameter a is 
given by Hertz' law F^'^^^ = B^^^"^, ^{t) = a — \fi — fj\ with B(a) being the 
elastic material parameter [13]. Assuming viscoelastic material properties the 
dissipative part of the contact force reads F^'^^^^ = A^/^^ [14-16] with the dis- 
sipative material constant A{a). Integration of Newton's equation of motion 
yields then [12, 17] 

e,{g) = 1 - D,g'/' + D,g'/' T . . . , (2) 



where the coefficients D1/2 depend on the elastic and dissipative particle ma- 
terial properties and on the radii of the particles (for details see [17]). 

The temperature of Granular Gases in the absence of spatial correlations is 
defined as the mean kinetic energy of the particles, just as for molecular gases. 
Due to dissipative collisions the temperature decays persistently and, there- 
fore, the thermal velocity vrit) = ^2T{t)/m decays too. For vanishing tem- 
perature the particles collide elastically, i.e., limT^o^ = 1- Since the formation 
of clusters in Granular Gases is a consequence of the dissipative particle colli- 
sions the question arises whether for gases of viscoelastic particles clusters may 
occur and persist. The linear stability analysis of the hydrodynamic equations 
of such gases supports the speculation that the cluster state is only a transient 
phenomenon for the case of viscoelastic particles [18]. 

For principal reasons (to be discussed below) the numerical investigation of 
the long time behavior of force free Granular Gases of viscoelastic particles 
is problematic. Therefore, in this paper we assume a coefficient of restitution 
with an extremely simplified impact velocity dependence: The particles collide 
with a constant coefficient of restitution e* if the impact velocity exceeds a 
certain value g*, otherwise they collide elastically: 



e{9) 



e* for q > q* 

^ ^ (3) 

1 ioi g < g* . 



We wish to mention that the coefficient of restitution of the form (3) was 
used in [19] as a numerical trick to avoid the inelastic collapse in event driven 



2 



MD simulations. In this context the transition velocity was chosen to be much 
smaller than the thermal velocity, g* ^ v^- Before investigating the cluster 
formation process first we will derive the law of temperature decay in the 
homogeneous cooling state. 



2 Decay of temperature in the homogeneous cooUng state 



In the absence of external forces a Granular Gas which is initialized at uniform 
number density n stays homogeneous during the first stage of its evolution, 
called homogeneous cooling state. In this regime the kinetic energy of the gas 
may be characterized by the granular temperature T", defined by the second 
moment of the velocity distribution function / {v) 

2^r=y ^/(^)c?^^, (4) 



(m is the mass of particles and d is the dimension) which decays persistently 
due to inelastic collisions. For e = const, the decay is described by Haff's law 
Tf/(t) [20], whereas for gases of viscoelastic particles the temperature is given 
by %{t) [17] 

TH(t) = o, TJt) = , (5) 



where th and t„ are the relaxation times [21]. 

Equations (5) have been first derived with the assumption of a Maxwell dis- 
tribution for the particle velocities. It is known that the velocity distribution 
function of Granular Gases deviates slightly from a Maxwell distribution [1- 
3], however, the functional form of Eqs. (5) keeps conserved also when we 
take into account these deviations. Before investigating the temperature de- 
cay let us first check whether the assumption of an approximative Maxwell 
distribution is justified for our case too. 

Starting at a certain temperature which corresponds to a thermal velocity 
vt S> g* almost all collisions occur with a constant coefficient of restitution 
£*, hence, in this range we expect to find in the homogeneous cooling state a 
velocity distribution close to a Maxwell distribution with the small deviations 
mentioned above. At late times when the thermal velocity is small with re- 
spect to the transition velocity, vt <^ g*, the majority of the collision occurs 
elastically. In this range we expect, therefore, that the particle velocities obey 
a Maxwell distribution just as for molecular gases. In the intermediate range 
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vt ^ g* Si sizable part of the collisions, i.e. collisions between slow particles, 
occurs elastically, whereas fast particles mainly undergo dissipative collisions. 
Therefore, it is not a priori clear whether the distribution function is close to 
the Maxwell distribution. 

We have simulated a two-dimensional force free cooling Granular Gas of = 
10^ particles of unit mass which collide according to the rule (3) by event- 
driven Molecular Dynamics. The parameters are: e* = 0.6, g* = 0.1, rj = 0.1, 
and the initial temperature is T{0) — 1 which corresponds to an initial thermal 
velocity vt{0) — V^. Figure 1 shows the scaled velocity distribution function 

m 

— * 



for the three mentioned cases together with the Maxwell distribution. Prom 
the results we conclude that even in the transition region vj- ^ g* the Maxwell 
distribution is a good approximation for the velocity distribution function. 




scaled velocity c 



Fig. 1. Scaled velocity distribution function obtained by MD simulations at times 
ti < t2 < ts where vrih) 3> g* (circles), VT{t2) ~ g* (squares), and vrits) ^ g* 
(plus symbols). The full line shows a Maxwell distribution. 



Using the standard approach based on the Boltzmann equation (e.g. [1-3, 
21]) with the assumption of a Maxwell distribution the equation for the time 
evolution of temperature may be obtained 

dT 2 

-^^-292{<yV-'nvTTf,2, (7) 



where a is the diameter of the particles, (72(0") is the contact value of the 
pair distribution function (5'2(c) = (1 ~ 77^/16) /(I — 1])'^ for a two-dimensional 
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gas with packing fraction rj = Tia^n/i), and fi2 is the second moment of the 
dimensionless coUision integral, which may be written in the form [3,21] 

IJ,2 = ^J dci J dc2 J deQ{-ci2 ■ e )/(ci)/(c2)(l - |ci2 • ef . (8) 

The integration is performed over the velocities Ci, C2 of the colliding pair and 
over the unit vector e = 6*12 as introduced in Eq. (1). The unit step function 
Q{x) guarantees that collisions occur only between approaching particles. 

Approximating the velocity distribution function by the Maxwell distribu- 
tion, /(c) = 7r~'^/^ exp(— c^), the above moment may be calculated. For two- 
dimensional systems it reads 



where u = u{t) = T(t)/To is the reduced temperature (Tq = T(0)) and go — 
9*/vt{0)- Introducing the variable x — 2u/gl we obtain from Eq. (7) 

i = -Q;(x^/2 + xV2)e-V- (10) 



where 



a 



^go{l-e*')92{a)anJ^. (11) 



If Vt{0) ^ g*, i.e. go 1, then a; ^ 1 at the initial stage of the temperature 
evolution. Equation (10) reduces then to i; = —ax^^'^, yielding Haff's law (5) 
with 



a 



90^ 2 



(12) 



For the late stage of the temperature evolution, i.e. for sufficiently small x, 
the general solution of Eq. (10) can be found as a series 



il + x) 



li'i + x) 1 (15 + 10,T + 3x2) 

1 -I -T A ^ 

2{l + x) 4 (l + x)2 



x^-r 



at + C, (13) 



where C is an integration constant. For large times, i.e., x <^ 1 it simplifies to 

x^/^eV^ = at (14) 
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Applying the (implicit) definition of Lambert's W function [22] 



(15) 



we find the analytical solution of Eq. (14) 
1 3 / 2 \ 2 T 



For large values of t we get the approximation (choosing the appropriate 
branch of W{x)) 



1 

X 



lng(at)^/3)+lnln(^(at)^/3) + ..; 



(17) 



which yields the asymptotics for t — > cxo 

T(t) = ^-^, (18) 
2 1ogQ;t' ^ ^ 



with the constants go and a defined above. 

At late times most of the particles collide elastically due to Eq. (3), there- 
fore, only collisions of particles whose velocities belong to the high-energy tail 
contribute to the decay of temperature. For this reason the asymptotic tem- 
perature relaxation is extremely slow as compared with the power laws Eq. 
(5), see Fig. 2. 




Fig. 2. Evolution of the granular temperature. Full line: MD simulation, dashed line: 
HafF's law (5), dash-dotted line: asymptotics due to Eq. (16). 
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For Granular Gases of particles which collide with e = const, it is known 
that the high energy tail of the scaled distribution function does not obey a 
Maxwell distribution (~exp(— c^)) but decays as exp(— const, x c) [4]. This 
overpopulation of the high energy tail arises due to a mechanism which is 
specific for gases of dissipative particles: the scaled velocity of a particle, c = 
v/vt can increase either by collisions with other particles or by keeping its 
velocity v while the thermal velocity vt, i.e. temperature, decreases. Both 
processes cause increase of the number of fast particles, i.e., they enhance 
the population of the high velocity tail. For very fast particles, c ^ 1, the 
former process may be neglected and the latter one leads to the exponential 
distribution function cxp(— const, x c) [4]. We wish to stress that for e ~ 
const, the Granular Gas has no inherent time scale and the shape of the 
scaled velocity distribution function /(c) does not depend on time. Hence, 
the exponential overpopulation of the high- velocity tail is a time-independent 
effect. 

For the coefficient of restitution according to Eq. (3) the overpopulation of 
the tail will not be observed at late times. The temperature decays only log- 
arithmically slow as compared with the power law for e — const. Therefore, 
the process of depopulation of the tail due to coUisions is not balanced any- 
more by the decrease of the temperature of the gas, which enhances the tail 
population. According to this argument in difference to the case e = const, we 
expect an underpopulation of the high-energy tail. There is a non-negligible 
probability that particles collide in such a way that the after-collision velocity 
of one of the collision partners exceed g* even if the pre-coUision velocities of 
both partners are smaller than g*. The probability of such collisions, which 
are responsible for the population of the velocity range v' > g*, decays rapidly 
as a function of the after-collision velocity v'. To estimate the affected range 
of particle velocities we consider the collision of elastic particles: The maximal 
after-collision velocity which may emerge in a collision of elastic particles with 
pre-coUision velocity g* is \/2g*. Hence we expect that only the interval of 
velocities {g*, V^g*) will be significantly populated. 

At present it remains unclear whether the deviation from the Maxwell dis- 
tribution at high velocities does change the behavior of the gas at late times 
significantly. This intriguing problem is subject of current research. 



3 Cluster formation 

The spontaneous formation of clusters in a force free cooling Granular Gas can 
be understood by simple arguments [8]: Consider fluctuations of the density 
in an otherwise homogeneous Granular Gas. In denser regions the particles 
collide more frequently than in more dilute regions, therefore, dense regions 
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of the system cool faster than dilute regions and the local pressure decays in 
these colder regions. The resulting pressure gradient causes a flux of particles 
into the denser region, which leads to further increase of the density. Hence, 
small fluctuations of the density are enhanced which leads to the formation of 
clusters. 

These arguments are certainly valid for the case e = const. For viscoelastic 
material properties, however, with decaying thermal velocity the collisions 
occur less and less dissipatively since approaches one with propagating 
time. Hence the question arises whether cluster formation in Granular Gases 
of viscoelastic particles occurs as a transient phenomenon only. This question 
is difficult and there is no conclusive answer yet, although the linear stability 
analysis [18] seems to support the hypothesis that cluster formation is only a 
transient phenomenon in Granular Gases of viscoelastic particles. 

The numerical simulation of chistcring Granular Gases is problematic since 
clusters grow rapidly and at a certain time their extension is comparable with 
the system size. From this instant on the simulation becomes invalid since the 
cluster interacts with itself via the periodic boundary conditions. For the case 
of viscoelastic particle properties e — e{g) we expect that clusters dissolve at 
late times when the typical coefficient of restitution e {vt) approaches one. 
To prove this by MD simulations we need a system which is large enough to 
assure that at no time the cluster size reaches the system dimensions. At the 
moment we are able to simulate systems up to about = 3 x 10^ particles, 
which seems to be far below the necessary size. 

We have simulated a system of n = 10^ particles which collide according to 
Eq. (3). This collision rule reveals the same asymptotic behavior as the rule 
(2) for the collision of viscoelastic spheres, lim^^o^^ = 1) nevertheless, the law 
is drastically simplified. Snapshots of the simulation are shown in Fig. 3. 

Starting with a uniform distribution after a short period of homogeneous cool- 
ing density inhomogeneities appear which grow up to a certain typical size. 
The collision frequency in dense regions is larger than in dilute regions which 
leads to decreasing thermal velocity. In the due of time more and more colli- 
sions inside the clusters occur with relative velocities g < g*, i.e., elastically. 
Therefore, the clusters dissolve. At late times the system returns to the ho- 
mogeneous cooling state. 

Why do we not observe system spanning clusters as in the case of £ = const.? 
The reason is the dynamics of cluster growth. To grow a cluster of a given size 
the system needs a certain time which is determined by the density and by the 
value of the coefficient of restitution. We start our simulation at a temperature 
which corresponds to a thermal velocity which is well above g*, i.e., the gas 
has essentially the properties of a Granular Gas with e — const. Therefore, 
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Fig. 3. Snapshots of a simulation of iV = 10^ particles (e* = 0.6, g* = 0.1, rj = 0.1). 
The pictures show the system after 0, 50, 100 (first row), 150, 200, 250 (second row), 
300, 500, 750 (third row), 1000, 1250 and 5000 (last row) collisions per particle. The 
initial temperature is Tq = 1 which corresponds to the thermal velocity vt = a/2- 
Starting with homogeneously distributed particles we observe cluster formation up 
to a certain cluster size. At late times the system returns to the homogeneous cooling 
state. 
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clusters begin to form and grow due to the described pressure instability. At a 
certain time the typical velocity inside the clusters falls below g*, i.e., from this 
moment on the clusters start to dissolve. This time corresponds to a typical 
cluster size which is significantly smaller than the system size. Therefore, we 
do not observe clusters beyond a certain typical size. 

According to a similar argumentation we expect clusters in Granular Gases of 
viscoelastic particles to dissolve too. Due to numerical limitations at present 
time it is not possible to perform the corresponding simulations without en- 
countering the problem of boundary conditions. 

In the recent study [19], where the restitution coefficient of the form (3) was 
used to avoid the inelastic collapse, the dissolution of clusters has not been 
detected. The reason for the persistent development of clusters in [19] was a 
very small value of g*, which was at least 100 times smaller than the thermal 
velocity of particles during the simulations. In our numerical experiments the 
clusters started to dissolve when the thermal velocity became comparable with 
g*. This regime has not been reached in simulations of [19]. 



4 Conclusion 

We have investigated a model Granular Gas of particles which collide with a 
stepwise coefficient of restitution which mimics the coefficient of restitution of 
viscoelastic particles: The particles collide elastically {s = 1) for small impact 
velocities and with e = const. < 1 if the relative particle velocity exceeds a 
certain value {g > g*). This simplified dependence of the restitution coefficient 
reflects the main feature of viscoelastic particles that coUisions tend to occur 
elastically for decreasing impact velocity. We analyze the temperature decay 
and the dynamics of cluster formation. We observe that clustering occurs in 
this model only as a transient process. The important property of the collision 
model is that the clusters dissolve before system spanning clusters appear. 
Hence, this model allows to investigate the total process without the drastic 
influence of the boundary conditions. The results of our study supports the 
hypothesis that clustering occurs as a transient process for realistic gases of 
viscoelastic particles, as was predicted before by means of the linear stability 
analysis. 



References 

[1] A. Goldshtein, M. Shapiro, Mechanics of collisional motion of granular 
materials. Part 1: General hydrodynamic equations, J. Fluid Mech. 282 (1995) 



10 



75. 



[2] T. P. C. van Noije, M. H. Ernst, Velocity distributions in homogeneous granular 
fluids: the free and the heated case, Granular Matter 1 (1998) 57. 

[3] N. V. Brilliantov, T. Poschel, Deviation from Maxwell distribution in granular 
gases with constant restitution coefficient, Phys. Rev. E 61 (2000) 2809. 

[4] S. E. Esipov, T. Poschel, The granular phase diagram, J. Stat. Phys. 86 (1997) 
1385. 

[5] N. V. Brilliantov, T. Poschel, Self-diffusion in granular gases, Phys. Rev. E 61 
(2000) 1716. 

[6] J. J. Brey, M. J. Ruiz-Montero, D. Cubero, R. Garcia-Rojo, Self-diffusion in 
freely evolving granular gases. Physics of Fluids 12 (2000) 876. 

[7] T. Poschel, S. Luding (Eds.), Granular Gases, Vol. 564 of Lecture Notes in 
Physics, Springer, Berlin, 2001. 

[8] I. Goldhirsch, G. Zanetti, Clustering instability in dissipative gases, Phys. Rev. 
Lett. 70 (1993) 1619. 

[9] R. Brito, M. H. Ernst, Extension of Half's cooling law in granular flows, 
Europhys. Lett. 43 (1998) 497. 

[10] F. G. Bridges, A. Hatzes, D. N. C. Lin, Structure, stability and evolution of 
Saturn's rings. Nature 309 (1984) 333. 

[11] T. Tanaka, T. Ishida, Y. Tsuji, Direct numerical simulatin of granular plug 

flow in a horizontal pipe: The case of cohesionless particles (in Japanese, for 
an english presentation of this work see [23]), Trans. Jap. Soc. Mech. Eng. 57 
(1991) 456. 

[12] R. Ramirez, T. Poschel, N. V. Brilliantov, T. Schwager, Coefficient of restitution 
for colliding viscoelastic spheres, Phys. Rev. E 60 (1999) 4465. 

[13] H. Hertz, Uber die Beriihrung fester elastischer Korper, J. f. reine u. angewandte 
Math. 92 (1882) 156. 

[14] N. V. Brilliantov, F. Spahn, J.-M. Hertzsch, T. Poschel, A model for collisions 
in granular gases, Phys. Rev. E 53 (1996) 5382. 

[15] W. A. M. Morgado, I. Oppenheim, Energy dissipation for quasielastic granular 
particle colHsions, Phys. Rev. E 55 (1997) 1940. 

[16] G. Kuwabara, K. Kono, Restitution coefficient in a collision between two 
spheres, Jpn. J. Appl. Phys. 26 (1987) 1230. 

[17] T. Schwager, T. Poschel, Coefficient of restitution of viscous particles and 
cooling rate of granular gases, Phys. Rev. E 57 (1998) 650. 

[18] N. V. Brilliantov, T. Poschel, Hydrodynamics of granular gases of viscoelastic 
particles, Phil. Trans. R. Soc. Lond. A 360 (2001) 415. 



11 



[19] X. Nie, E. Ben-Naim, S. Y. Chen, Dynamics of freely cooling granular gases, 
Phys. Rev. Lett. 89 (2002) 204301. 

[20] P. K. HafF, Grain flow as a fluid-mechanical phenomenon, J. Fluid Mech. 134 
(1983) 401. 

[21] N. V. Brilliantov, T. Poschel, Velocity distribution of granular gases of 
viscoelastic particles, Phys. Rev. E 61 (2000) 5573. 

[22] R. M. Corless, G. H. Gonnct, D. E. G. Hare, D. J. Jeffrey, D. E. Knuth, On the 
Lambert W Function, Adv. Comput. Math. 5 (1996) 329. 

[23] Y. Taguchi, Powder turbulence: Direct onset of turbulent flow, J. Physique 2 
(1992) 2103. 



12 



